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Abstract 

We discuss the features of a two-Higgs doublet model exhibiting a two stage 
phase transition. At finite temperatures electric charge violating stationary points 
are developed. In conjunction with CP violation in the Higgs or the Yukawa sector, 
the phase transition to the charge conserving vacuum, generates a net charge 
asymmetry AQ, in the presence of heavy leptons, which may be well above the 
astrophysical bounds put on AQ unless the heavy leptons are sufficiently massive. 
This type of transition may be of relevance for supersymmetric extensions of the 
Standard Model, since it shares the same features, namely two Higgs doublets and 
similar CP violating sources. 
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A major problem in cosmology is the explanation of the observed baryon asymme- 
try One of the great achievements of the standard hot big bang cosmological model 
is certainly nucleosynthesis, which explains the measured abundances of the light el- 
ements in the Universe. It leads, however, to a fine tuning problem of one parame- 
ter, the baryon to entropy ratio. It is required that this parameter be in the range 
4 x 10~ n < ub/s < 7 x 10~ n , where % is the net baryon density in a commoving 
volume and s is the entropy density. 

The requirement for cosmological baryon asymmetry poses severe constraints on the 
allowed particle models. In order to generate the baryon asymmetry dynamically, as the 
Universe cools, the following criteria should be satisfied: The baryon number is not 
conserved, C and CP symmetries are violated and there are non- equilibrium conditions. 

The only successful attempts so far, in explaining baryogenesis at the electroweak 
scale, are those based on the extensions of the Standard Model (SM) with two-Higgs 
doublet models. The CP violation provided by the Kobayashi-Maskawa matrix proves 
to be too small to explain cosmic baryon asymmetry 0. These models have a new 
source of CP violation, the phase between the two vacuum expectation values of the 
Higgs fields Shaposhnikov |7J and MacLerran H first proposed that the two-Higgs 

model is a candidate model to explain electroweak baryogenesis. Turok and Zadrozny 
subsequently analyzed how it could work |§. 

In supersymmetric models it is necessary to have at least two Higgs doublets, while 
extra CP violating sources, occurring in other sectors, arise naturally. In this case the 
magnitude of CP violation is severely constrained by neutron's dipole moments and/or 
unification assumptions, restricting the allowed parameter space [|T0|^T^1 . 

In the two-Higgs model the desired baryogenesis at the electroweak scale is imple- 
mented by the introduction of terms that break CP invariance explicitly at the tree 
level. The necessity for introducing such terms is the natural suppression || of Higgs 
mediated flavour changing neutral currents (FCNC). In fact, it is well known, that CP 
symmetry is spontaneously broken in the presence of two Higgs doublets, at the expense 
of introducing unnatural suppression of FCNC processes \TE\. On the other hand, CP 
violation with two Higgs doublets and FCNC processes that are suppressed due to the 
heaviness of the Higgses involved [TO , lead to second order phase transitions making the 
baryogenesis scenario almost impossible to implement [|T7|] . 

In the context of a two-Higgs model and in the absence of extra CP violating terms, 
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domains of the Universe can first tunnel to a new type of minimum in which either of the 
two Higgs doublets develop no vacuum expectation value, before evolving to the usual 
electroweak symmetry breaking minimumF]. In other words the phase transition occurs 
in two stages. The first stage has been shown to be a first order transition, while the 
subsequent is of second order fl9[]. Fluctuations of the minimum, which we have due 
to the presence of the extra terms, are then proved to lead to amplified CP violation 
at high temperatures ]19| in the same spirit with Ref. |2^,|2l|, resulting to reasonably 
large baryon asymmetry. This happens in a significant range of the parameters that may 
cover the Minimal Supersymmetric Standard Model (MSSM). 

The present work is an analysis of the two-Higgs phase transition in the sense that we 
fully take into account the contribution of the extra breaking terms and do not consider 
them as small perturbations. The calculation reveals a very interesting feature: The 
phase transition takes place again in two stages. This happens in the entire range of 
the parameters, as long as the ratio of the vev's is not fine tuned to unity. The explicit 
CP violating angles developed in this first stage, are amplified compared to the zero 
temperature ones, as required for baryogenesis. During this stage a stationary point 
that breaks the U(l) gauge symmetry of electromagnetism is always developed. In the 
context of the two-Higgs extension of the SM the charge violating cosmological phase 
does not produce a net charge asymmetry. However if the spectrum contains heavy 
leptons, as we will discuss, it is possible to create a net charge asymmetry. Such leptons, 
could be heavy Majorana neutrinos. A strong motivation of having heavy neutrinos is 



to explain the smallness of neutrino masses using a see-saw mechanism . 

The presence of a charge asymmetry can lead to dramatic effects in the cosmological 
evolution due to the large strength of the electromagnetic interactions as compared 
to those of gravity. Any specific model that produces a charge asymmetry beyond a 
certain limit is ruled out on cosmological grounds. Thus only models yielding small 
charge asymmetries are acceptable. It should be noted that small values for the charge 



asymmetry are able to produce a primordial magnetic seed field [2t]. 

Using weak isospin doublets both of weak hyperchargef] Y wcak = +1 the two-Higgs 



1 Land and Carlson ]18| first noticed the possibility of a two stage phase transition. However the 
characteristics of it are different from this we are discussing here. 

2 We follow here the notation of Ref. p7J in which both Higgs doublet fields have same hypercharge. 
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scalar potential can be written as follows |27| 



^Hig gs = ^l^ + ^^ + A^lf^ + A^^^+Aalfi^)^) 

+A 4 ($ t 1 $ 2 )($^ 1 ) + ^[(S^) 2 + (4®i) 2 ) + Vb , (1) 

where Aj are real numbers and 

$ 1 = ^f 01+ ^ 2 V $ 2 = * ^ 5 + ^V (2) 
y/2\(h + i<l>J y/2\<h + i<l>J 

The above potential, with the exception of Vd which we discuss in the following, is 

the most general one satisfying the following discrete symmetries: 

$ 2 -> -$ 2 , $i -> $i, <f fl -> -<f fl , It^j -> < , (3) 

where u l R and represent the right-handed weak eigenstates with charges | and — | 
respectively. All other fields involved remain intact under the above discrete symmetry. 
This symmetry forces all the quarks of a given charge to interact with only one doublet, 
and thus Higgs mediated flavour changing neutral currents are absent. When this discrete 
symmetry is broken, during a cosmological phase transition, it produces stable domain 
walls via the Kibble mechanism ||28|| . One can overcome this problem by adding terms, 
which break this symmetry providing at the same time the required explicit CP violation 
for baryogenesis. The most general form of that part of the potential which breaks this 
discrete symmetry is 

Vt> = -yU3$I$2 + A 6 ($ t 1 $ 1 )($ t 1 $ 2 ) + A 7 (^$ 2 )($ t 1 $ 2 ) + h.c. (4) 

This will shall call hereafter as D— breaking part. The parameters fa, Xq and A7 are in 
general complex numbers 

= ml e ie \ A 6 = k e id «, A 7 = l 7 e i(>7 (5) 

providing explicit CP violation at the tree level. 

In order to study the structure of the vacua we can perform an SU(2) rotation, 
that sets the vev's of the fields for 01,2,4 of Eq. (H) equal to zero. Solving the system 
dV/dipi = implies several different stationary points. One of them is the usual asym- 
metric minimum, that respects the U(l) of electromagnetism 

*, * (°V * * ( V ( 6 ) 



In Eq. u, v, <p are real numbers. The phase tp is the explicit CP violating angle that 
appears due to existence of the D— breaking termsQ. The acceptable parameters of the 
model are those ensuring that the above stationary point becomes the absolute minimum 
at zero temperature. 

The free parameters of the model can be taken to be the quartic couplings Aj, the 
ratio (3 = u/v and the mass parameter fis, since we have the following conditions at the 
stationary points 



f4 
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(A 3 + A 4 + A 5 ) v 2 + ml cos 9 3 f3 1 — m\ sin 9 3 (3 1 (p 



11 3 3 

-~l 7 cos 9 7 v 2 fi~ x + ~l 7 sin 9 7 v 2 (3~ 1 ip — -l e cos 9 e v 2 (3 + —l 6 sin 9 e v 2 /3(p, (7) 



and 



l4 = -X 2 v 2 - - (A 3 + A 4 + A 5 )/?V + m^cos^ 3 P 

3 1 1 

-~l 7 cos 8 7 v 2 f3 + -l 7 sin 67 v 2 (3cp - -l 6 cos 9 6 v 2 f3 3 (8) 

4t Zi 4r 

— m\ sin 63 + jl 7 sin 6 7 v 2 + sin 9 6 v 2 (3 2 
X5V 2 [3 — m\ cos 6*3 + \l 7 cos 6 7 v 2 + -Aq cos 6q v 2 (5 2 

In deriving Eqs. (0) and @ we have made use of the fact that the zero temperature 
phase is small \ip\ 1, as discussed previously. The value of ip can become very small 
either by taking the phases 9 3 ,9q, 9 7 to be small or by choosing the quartic coupling A5 
to be large enough. One should choose the free parameters ensuring that the stationary 
point of Eq. (^) is indeed a minimum and the potential is bounded from below. 

Scanning the whole parameter space would be too time consuming. Therefore in our 
analysis we consider values of the quartic coupling constants such that the ratios 
with % = 1,2, ...,5, vary in the range from 0.01 to 1.00, where g is the SU(2)l coupling 
constant. The set of the parameters displayed in Table [l] is a sample which respects the 
experimental limits on the Higgs masses, the condition that the potential is bounded 
from below, and that there is an absolute U(l) em minimum at zero temperature. The 
couplings A 6 and A 7 should be taken sufficiently large to respect the experimental limits 
from Higgs searches. 

3 The value of the phase ip is severely constrained by data on neutron's dipole moment and will be 
taken small in the following. 
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The cosmological phase transition of the model under consideration can be studied 
using the finite temperature effective potential. A handy feature of two-Higgs doublet 
models is that they can give a stronger first order phase transition and thinner bubble 



walls 29-31 



In order to get a complete picture of the transition it is necessary to explore the 
full range of the potential. At T = the right vacuum is chosen to be the absolute 
minimum after symmetry breaking takes place. For a model with two complex doublets 
this is realized by selecting the two neutral components to develop non-vanishing vacuum 
expectation values. In other words the minimum rests in a specified plane. However this 
by no means ensures that during the phase transition the absolute minimum remains 
in the same plane, since the field space is effectively five-dimensional. What is usually 
assumed is that the minimum remains in the same plane, as in the zero temperature 
case. However this picture may not be correct, leading to misleading conclusions. 

In our study we include the one-loop radiative corrections to the tree level potential, 
using the temperature corrected fermion and gauge boson masses^. Expanding the loop 
integrals we get a cubic in the scalar fields term, which stems from the Matsubara 
zero mode |]32j]. It is beyond the scope of this paper to discuss well known problems 



associated with the validity of the perturbation and the high temperature expansion. 
Thus we employ the simplified potential 

V(T) = V mggs + I E(Mi), + 2m?]T 2 - i- ^ M a)T T, (10) 

i i 

8 

where ml — <fij, is the field dependent top quark mass and (M|)j are the corre- 

i=5 

sponding eigenvalues of the gauge boson mass matrix 

2 

(M\) ab = g 2 ^& k T a T b <& k . (11) 

k=l 

The Lie algebra matrices appeared in Eq. (|ll]) are defined as follows 

T a = a a for a = 1, 2, 3 , 

T a = tl for a = 4, (12) 



4 For simplicity wc ignore the contribution of the scalar bosons. Their presence makes the linear in 
temperature term smaller and the first order transition weaker. However the main results of the paper 
regarding the charge breaking intermediate phase, remain the same irrespectively of the specific form 
of the effective potential. 
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with t = g'/g and g (g r ) are the SU{2)l {U{1)y) coupling constants. 

In order to find the location of the stationary points of the potential we adopt the 
following procedure: At any given temperature and for given set of couplings we first 
check the shape of the potential in every two dimensional plane (0j, <pj) which is specified 
by <pk = 0, k 7^ i, j, in order to locate all possible stationary points. These can be either 
extrema or saddle points. Their nature is identified by simply looking at the signs of the 
eigenvalues of the second derivative of the potential. 

What the analysis reveals is stated in the following: 

• When (3 = 1 and at very high temperatures the symmetry is unbroken, as expected. 
As the temperature drops, the first asymmetric minimum develops. This is charge 
preserving, see Eq. (H), and lies on the plane 03 = 07, 08 7^ and 0j = 0. There 
is a barrier between this minimum and the symmetric one, which below a critical 
temperature lies higher, so we have a first order phase transitionF]. In this case 
everything looks conventional. However we should point out that the situation of 
having (3 = 1 is unnatural in the sense that a fine tuning of the parameters is 
required in order for this to be realized. 

• For (3 < lthe picture alters considerably^. In these for instance the 
one displayed in Table [l], an additional stationary point which is charge breaking 
appears. Depending on the values of the parameters involved it may have appeared 
earlier or later than the charge conserving one. For a qualitative discussion of the 
thermal evolution of the Universe it suffices to distinguish three subcases: 

(i) At very high temperatures only the symmetric vacuum Vs appears. As temperature 

drops another local minimum Vex-, which is charge violating, starts developing 

which below a critical temperature lies lower than the symmetric minimum (Vex < 

Vs). During this stage bubbles of this minimum nucleate and propagate. As the 

Universe further cools a stationary point, Vc, starts formatting which is charge 

preserving but lies higher than Vex- At this stage this is a saddle point, but as the 

temperature further drops by about AT ~ 100 GeV, Vc moves lower and becomes 

the absolute minimum. At the same time Vex becomes unstable and through a 

second order phase transition the scalar fields roll towards the stable minimum Vc- 

5 However this does not guarantee bubble formation with kink- like walls |33[) . 
6 The case (3 > 1 is identical to (3 < 1 with the interchange $i ^ $2- 
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The kind of evolution just described occurs for instance when (3 < 0.8, A6,7 > g 2 , 
m 3 100 GeV but the bound on can be relaxed if the parameter (3 is very 
small. However this is not the only parameter region where this interesting subcase 
is realized. 

(ii) There are also values in the parameter space for which, after the symmetric vacuum, 
the charge preserving minimum Vc is developed. At a later stage of the evolution, 
a charge violating saddle point Vex starts forming which stays higher than Vc 
during the whole evolution. Then the final stage of the evolution in this case does 
not differ from the subcase discussed previously. A region of the parameter space 
where this holds is for 0.8 < (3 < 1, Xqj > g 2 , m 3 > 110 GeV, but as in the 
previous case, we can find other regions too where this subcase is realized. 

(iii) There is also a parameter range where both stationary points exist (of course one 
of them is the minimum and the other is a saddle point) with comparable depths 
and comparable barriers from the symmetric vacuum. Therefore transition to these 
points proceeds with similar tunneling rates leading to the nucleation of both types 
of bubbles for a while. It is worth mentioning here, that the tunneling rate for the 
transition from the symmetric vacuum to the aforementioned stationary points 
receives contributions from paths not belonging to the plane of the two points, 
due to the fact that the potential is multidimensional. These corrections, which 



are given in a form of a determinant |34|] in front of the well known exponential 



law |[35|| , although seem to affect little, there are cases where their contribution 
is significantly enhanced p6| . In the simple cases, it suffices to use the one-loop 



corrected tunneling rates [37 



In order to study the minima of the potential at finite temperature in the eight 
dimensional field space, we exploit the SU(2) invariance and eliminate three out of the 
eight degrees of freedom. In fact by an SU(2) rotation R(T) = e" 7 "-^^ the doublets 
$i, $2 can be brought into the following standard form: 



where u(T) is real and z(T), w(T) complex. The reduction to five component simplifies 
the picture a great deal. With the above parameterization the two types of the stationary 
points discussed previously, are as follows: 
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Charge conserving: u(T) ^ 0, z(T) = 0, w(T) = v(T)e lv( - T * ) ^ (where v is real) , 

Charge breaking: u{T) — 0, z(T) ^ 0, w(T) = v'{T)e iv '^ ^ (both 2, u> complex) . 

There is a degeneracy in the vacuum structure of the charge breaking minimum, since 
the potential V(T) in Eq. ( [T0| ) has an SU(2) symmetry, if one sets $1 = 0. Therefore, 
there is a manifold of absolute minima satisfying 0j = 0, for i — 1, 4 , and 4>\+4>1+<^+ 
02 = \ Z (T)\ Z + \w{T)\\ This vacuum degeneracy forbids the creation of net baryon and 
charge asymmetry p0| , |21|j . What we should stress at this point is the significance of the 
D— breaking terms. Their presence signals the existence of charge breaking stationary 
points. 

Another thing that the multidimensional analysis of the phase transition makes clear, 
is that the explicit CP violating phase <p{T), at high temperatures is significantly en- 
hanced for all the examined sets of parameters, as compared to its zero temperature value. 
Recall that <f{T) is non- vanishing, provided that the D— violating terms of Eq. (f|) are 
present. This is true even in the case where the phases $3,6,7 are quite small. This is 
due to large temperature depended factors in front of the arguments of the D— breaking 
terms |L9| . 



The development of a large <p(T) affects the baryogenesis scenario, since the baryon 
asymmetry depends linearly on (p(T). In fact within the context of the popular non-local 
mechanisms, which seem to work more efficiently, the fermion which is reflected from 
the bubble wall experiences a space dependent phase. Since right-handed fermions have 
different reflection coefficients from the left-handed anti-particles a baryon asymmetry 



is produced P3[, which for small velocities v w of the bubble wall is given by 



-v w f — )<p , (14) 



s 2g s 7i A \T C 

where / is a Yukawa coupling and g s denotes the number of spin states. 

A net excess of charge density pq = e Uq, is produced during the first charge breaking 
stage of the phase transition, if Sakharov's last two conditions hold as in the baryon case. 
If C and CP are conserved, then C and CP conjugate reactions produce equal amounts 
of opposite sign charge asymmetry resulting to zero net charge asymmetry. As we shall 
see, in our case we have CP violation in charge violating processes, resulting either 
from the D— breaking terms in Eq. (j4|), which are in general complex, or from complex 
Yukawa couplings. The number density of particles and antiparticles is zero in thermal 



S 



equilibrium, since the masses of a particle and its CPT conjugate antiparticle are equal. 
However this is not the case during a phase transition. 

The idea of a charged Universe was first considered by Lyttleton and Bondi [[^J . The 
fact that the strength of electromagnetic interaction is much larger that the gravitational 
one, causes any net charge asymmetry to have serious cosmological implications. Soon, 
upper limits on the magnitude of the allowed asymmetry were posed (for recent results 



see Refs. |40|, f4l| ) . The most conservative constraint comes from the requirement, that 
the gravitational interaction is larger than the electrostatic repulsion of the net charge. 
This gives uq/s < 10 -18 , where uq is the number of elementary charges in excess. The 
entropy density is used as a measure of the commoving volume, since from entropy 
conservation we get s oc R~ 3 . 

More stringent bounds come from the demand that the anisotropic angular distri- 
bution of cosmic rays, induced by the electric field of the net charge, be lower than the 
observed one. This leads to ng/s < 10~ 30 . 

An additional bound is furnished by the observed isotropy in the cosmic microwave 
background radiation. A small charge asymmetry can cause very large anisotropies [|IJ , 
through the same mechanism of the gravitational density perturbations. Thus, charge 
asymmetries induce temperature fluctuations in the microwave background (Sachs- Wolfe 
effect). In this case the constraint is: uq/s < 10~ 29 . Any specific two-Higgs doublet 
model which does not respect the aforementioned bounds has to be excluded on cosmo- 
logical grounds. 

Now we will shall see that within the context of the SM, no net charge asymmetry 
can be produced, even if the Universe passes from a charge violating phase. However 
the case is very different in the presence of heavy leptons. 

As we have already discussed the phase transition happens in two stages. The U (l) em 
gauge symmetry can be broken during the first one. One can realize the violation of the 
U(l) em gauge symmetry because of the development of non-vanishing vev's in the charged 
components of the Higgs doublets. The Higgs part of Lagrangian of the two-Higgs SM 
model is: 

£ H i ggs = Yl (D^iYiD^) + £ Yuk - V mggs , (15) 

i=l,2 

where the Yukawa terms for one family of fermions are 

£vuk = -h u Q$ 2 u R - h d Q$id R - h e L$ x e R + h.c. , (16) 
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and Vmggs reads from Eq. ([!]) . Q and L denote quark and lepton SU (2) l doublets, 
respectively. The Higgs doublets in Eqs. ( |T5| ) and (|TB[ ) are as follows 

^ = ( il 2 = ( 1- 1 • (17) 



1,2 / \ 



During the charge breaking phase only the second Higgs doublet $2 develops non- 
vanishing vev's, as have been discussed previously. The development of charged vev's 
has two major implications: (a) There is a mixing between "neutral" and "charged" 
gauge eigenstates resulting to new mass eigenstates without electric charge assignment, 
and (b) there are "electric charge" violating interactions among the gauge eigenstates. 

The presence of the charged vev's alters significantly the mixing of gauge bosons 
W ± ,W < ' 3 \ B []. In place of the usual mixing between the neutral gauge eigenstates 
W^ 3 \B, which results to the mass eigenstates of j,Z, there is a mixing between the 
neutral and charged gauge eigenstates (W^ 3 \ B, W /± ), yielding new mass eigenstates 
Vi, % — 1, ...,4. These eigenstates do not have an electric charge assignment. The second 
novel feature is the existence of charge violating interactions among gauge and Higgs 
bosons. 

The Yukawa part of the Lagrangian, after shifting Higgses from their vev's, involves 
the following bilinear terms: 

Auk = -Ku L UR{(j)lY - h d d L d R {(j>\) - Ke L e R {4>\) 

+h u d L u R {(f>2) - h d u L d R ((p+) - h e v L e R ((f)^) + h.c. (18) 

The first three terms yield the usual mass terms for the fermions, while the rest lead to a 
mixing between fermions of different electric charge. There are no additional interactions 
in the Yukawa sector, except those appearing in Eq. (|T8|). 

In the following for definiteness we shall elaborate the (3 < 1 case. Analogous results 
can be obtained for the (3 > 1 case. From Eq. (|T8| ) it is seen that only d and u quarks 
suffer a mixing. The mixing between e and v is absent, since (4>x' ) = 0. In this case the 
gauge and mass eigenstates of leptons are the same. 

As far as the potential is concerned there are new mass terms and interactions, that 

mix charge and neutral components of Higgs doublets. These are of two types: (a) terms 

7 The are the usual linear combinations -j= {W^ 1 ' =F iW^). During the charge breaking stage 
the ± labels do not correspond to any physical electric charge assignment. 
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that arise from the D— violating part of the potential, which carry CP violating sources 
(A*3, ^6,7)) an d (b) terms that arise from the rest quartic part, which are CP preserving. 
Terms carrying CP violating phase will be proved extremely important in producing net 
electric charge during the first stage of the phase transition. 

We intend to calculate the net electric charge production, during the U(l) em breaking 
stage of the phase transition. Therefore we focus our attention on CP-odd charge vio- 
lating interactions, that participate in charge violating processes. CP violating sources 
are the complex phases in the Higgs potential (D— violating terms couplings) and in the 
Yukawa terms (Yukawa couplings). Although the number density of the light leptons is 
abundant in the plasma, leptonic processes, such as for instance the decay I —>■ v\ vv, 
yield vanishing net electric charge, because there is always sufficient thermal energy to 
activate the inverse reaction, v\ vv — > I . Therefore the two processes coexist in equi- 
librium in the plasma, resulting to a vanishing net electric charge. Thus we turn our 
attention on these charge violating processes, in which the initial state is characterized 
by a mass which is large, as compared to the masses of the particles produced in the 
final state. In these cases only when the produced particles acquire large kinetic energies, 
due to their thermal motion, the inverse reaction can be activated. On these grounds a 
candidate process, within SM, to succeed in producing non-zero net charge asymmetry, 
is the charge violating decay t —>■ d(s) vv. However even in this case, the net effect is 
zero since there are fast strong reactions tg ^ d(s) g, which keep in equilibrium the 
quarks, washing out any asymmetry. Thus within the SM it is not possible to gain any 
net charge asymmetry. However in the case we have heavy additional particles which are 
out of equilibrium and participate in charge violating reactions, a net charge asymmetry 
can be produced. Such candidate particles are heavy leptons. 

As an example of how a charge asymmetry may be produced we consider a simple 



and favorable model f42H45f, in which the fermionic spectrum of the SM is enlarged by 
adding two isosinglet neutrino fields, n c and S, per family /. Their mass mixing, in Weyl 
basis, is given by 

1 (° D \ (" \ 

- Cmass = - ( v n c S ) \ D T M T rf + h.c. (19) 

1 \ M / \ S J 

The fields v, n c , S represent collectively all families, v\ = ( v e v^ v T ) so that D and 
M are actually 3x3 mass matrices. We assume that the elements of D are much smaller 
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than those of M. This model leads |^2| to three massless neutrinos u iL , i = 1,2,3, one 
for each family and three Dirac heavy neutrinos N a , a = 5,6,7 with masses ~ 
M + 0(D 2 /M). The mixing can be expressed as follows 

V lL = ( K Lhv'iL + Yl ( K H)laN aL , (20) 

i=l,2,3 a=4,5,6 

with Kh ~ 0{jj) and K L ~ 0(1). We can also introduce a small Majorana mass /i in 
the (3,3) entry of the mass matrix if we want the neutrinos v\ to acquire a small mass. 

The relevant interactions able to produce charge asymmetry are the Yukawa and 
gauge interactions 

C Y = - {h N f ( VZ , T L ) i ^ 2 (n R ) . + h.c. , (21) 

where (/iat) 1 "' are Yukawa couplings. The Si field does not couple to $j. The charged 
current part of the gauge interactions is given by 



2^2^ ^ 



(1 - 75) (K L ) U v[ + J2h, (1 - 75) (K H ) la N a 

a 



+ h.c. 



(22) 

Then we consider the process iV — > l<pi 2 , presented in Fig. [I], and for simplicity 
we assume one family. In order to calculate the net electric charge production rate we 
compare, as in the case of baryon asymmetry [ 17], |6] , the decay rate of process N — > I <p\ 2 
and that of its CP conjugate. The difference of these rates will be proportional to 
the produced average net electric charge. In Figs. filled boxes (circled crosses) are 
used to indicate CP-even (odd) charge violating mixings or interactions between the 
gauge eigenstates. The circled crosses depend on the complex couplings /x 3 , A 6; 7. All 
lines in Figs. |],|] represent gauge eigenstates, just in order to make apparent the charge 
violating mixings and interactions. These gauge eigenstates are linear combinations of 
the corresponding mass eigenstates. 

In Fig. 0a the tree level graph does depend on the CP violating complex parameter 
A7. In the Born approximation there is no net electric charge production, since the cross 
section remains invariant under the CP transformation, even if complex couplings are 
involved: 

r Bom (iv -> 1 0°) = r Born (]v -> . (23) 

The one- loop diagram displayed in Fig. |l|b depends linearly on the Yukawa coupling hi, 
which can be complex in general, and on A3 which is however real. The complexity of A7 
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and Yukawa couplings results to non-vanishing net electric charge production from the 
interference term of the tree level graph in Fig. la and the one-loop graph in Fig. [T|b. 
The contribution of this interference term can be written as: 



-unterf 



(N I 0?) = \ 3 h t I Nl | (fo) | X 7 (h N ) 2 + h.c. , (24) 



where Ini is the relevant kinematic factor corresponding to this interference term. For 
simplicity we have assumed that is real. The CP dependence of these elements is 
of the order of the phases of hi and A7. The CP conjugate process N — > 1 yields an 
interference contribution: 



-dnterf 



(N -> UT) = A 3 ^r I*m |(02 > I KM 2 + h.c. , (25) 



where in both Eqs. (gjj) and (|25|), we have used the modulus of (02 ) , since spontaneous 
CP violating phases cancel each other. 

The net electric charge produced during the first stage of the phase transition through 
the process N — > I (f)\ is proportional, at the one-loop order, to the quantity 

r interf (iv ^ /0 O ) _ r intc rf(7? ^ l0 O* ) ^ 4A3 |^| I m [ ( 26 ) 

To leading order, in the small CP violating phases, the following approximation holds 

Im[M 7 ] ~|A 7 |N(# 7 + arg(/iO) . (27) 

The average net electric charge produced via the above charge violating process, mea- 
sured in units of the electron charge, is therefore 

AQ = — [r intorf (N -»• 1 0?) - r interf (F -> 70°*)] , (28) 
where T^r ~ -^-M. Using Eq. fl26|) , Eq. fl28| ) can take the form 

AQ^^^ |<0 2 >||A 7 | N(0 7 + arg(fc,)) Im [J M ] • (29) 

For the process iV — » Z0?,, presented in Fig. [l|c-[l|d, we get a similar contribution to AQ 
with A3 replaced by A 2 . So the net charge density produced during the charge violat- 
ing phase, if all the CP violating phases are of order 10~ 4 as required by baryogenesis 
scenarios, turn out to be 

1(T 8 

AQ ~ — j- , (30) 

Mj TeV V ; 
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leading to 

AQ n N f 1(T 8 \ n N 
n Q/ s ~ —TZ— ~ ^/rrw TTT- ( 31 ) 



<7* n 7 \ M/ Te V / g* n 7 
For M = 1 TeV, ~ 10~ 5 , yielding nq/s ~ 10~ 13 which is clearly much larger than 
the bounds 10 _29 -10 -30 quoted previously. However the heavy neutrino density n^s falls 
rapidly with increasing M resulting to a smaller charge asymmetry. In order to find 
r^N at the electroweak scale T = 100 GeV one has to solve the appropriate Boltzmann 
equation. By solving this, we find that the bound 10~ 30 on hq/ s is saturated for a mass 
M ~ 5 TeV. For larger values of M the aforementioned upper bound on nq/s is always 
satisfied. Thus in such models the cosmological upper limits on the charge asymmetry 
impose lower bounds on the heavy neutrino masses, unless the combination of phases 
appearing in Eq. ( ^9|) is fine tuned to values much smaller than 10~ 4 . 

There are also other charge violating processes N — > I V (V denotes the neutral gauge 
bosons) that can produce charge asymmetry through the interference of the tree-level 
graph, presented in Fig. |2|, with the corresponding one-loop graphs. This process is less 
significant than the one we considered above since its rate is ^f- times smaller than the 
rate of the previous interaction. However if 9 7 = arg(/^) = then the process N ^ IV 
is the dominant mechanism for charge asymmetry provided 8 3 , # 6 are different from 
zero. When the the Yukawa coupling Hn of the heavy leptons is complex, the process 
N -> / cf>1 2 can lead to non-zero charge asymmetry if we adopt the model of Ref. j|5| 



and use a similar interference pattern with that presented in Ref. [f47 |. 

To summarize, in this paper we discuss the phase transition of a model with two 
Higgs doublets, in which CP violating terms are present in the scalar potential. The 
motivation for this study originates from the fact that the successful mechanisms for 
explaining baryogenesis are based on such models. Besides this, such an extended Higgs 
sector resembles that of supersymmetric (SUSY) models and hence our conclusions may 
be of relevance for SUSY extensions of the SM. 

Using the finite temperature one-loop corrected effective potential we found that the 
phase transition occurs in two stages. During the first stage of the transition CP violating 
angles are amplified, as required by baryogenesis, and simultaneously a charge breaking 
stationary point is developed which is a minimum, in a wide range of the parameter space. 
This leads to a non- vanishing charge asymmetry, in the presence of heavy leptons, due 
to the appearance of CP violating sources within the D— breaking terms and/or Yukawa 
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couplings. In the context of the SM no net charge asymmetry is produced. However in 
models with heavy leptons the magnitude of the asymmetry is found to exceed existing 
astrophysical bounds, constraining the mass spectrum of heavy leptons. In particular 
in models with heavy neutrinos we find that heavy neutrino masses smaller than 5 TeV 
result to unacceptably large charge asymmetry. 

The study undertaken in this paper maybe of relevance for supersymmetric models 
and other extensions of the SM which are characterized by a complicated Higgs sector 
and CP violating sources in the SUSY breaking terms of the scalar potential. The de- 
velopment of a charge breaking minimum during the cooling down of the Universe may 
lead, in this case too, to a large net charge asymmetry restricting the allowed parameter 
space. Besides this, one has to pay special attention to the appearance of other kinds 
of minima, such as color breaking, which develop at finite temperatures. These may im- 
pose further constraints on supersymmetry parameters and affect the phenomenological 
predictions. Such a study is under consideration and the results will appear elsewhere. 

It is worth mentioning that the possible two stage transition, studied in this paper, 
may be a feature shared by other multiscalar potentials deserving special attention. In 
GUT theories the breaking of U(l) em symmetry at the unification scale offers a working 
mechanism for the resolution of the monopole problem as proposed by Langacker and 

Pi ra. 
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Table 1: Sample of allowed parameters respecting all theoretical constraints and exper- 
imental bounds on Higgs boson masses (see main text). 
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Figure 1: Diagrams which contribute, up to one- loop order, to charge violating 
process N — > l(f>i 2 - Filled boxes (circled crosses) are used to indicate CP-even 
(odd) charge violating mixings or interactions between the gauge eigenstates. 




Figure 2: Tree level diagrams which contribute to charge violating process 
N — > I V. The filled box is as in Fig. |l|. 
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